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Abstract 

Using 0.8 x 10 6 D + D~ pairs collected with the CLEO-c detector at the ^(3770) resonance, 
we have searched for flavor-changing neutral current and lepton-number-violating decays of D + 
mesons to final states with dielectrons. We find no indication of either, obtaining 90% confidence 
level upper limits of B{D+ -> n+e+e-) < 7.4 x 10~ 6 , B{D + -► Ti~e + e + ) < 3.6 x 10~ 6 , B{D + -» 
K+e+e-) < 6.2 x 10~ 6 , and B(D+ -» R-e + e + ) < 4.5 x lO" 6 . 
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Searches for rare-decay processes have played an important role in the development of 
the Standard Model (SM). The absence of flavor-changing neutral currents (FCNCs) in kaon 
decays led to the prediction of the charm quark |l| , and the observation of B°-B° mixing, a 
FCNC process, signaled the very large top-quark mass 0. To date, rare and forbidden charm 
decays have been less informative and less extensively studied. In this Letter we present 
searches for the FCNC decays |EJI| D + — > n + e + e~ and D + — > K + e + e~, and the lepton- 
number-violating (LNV) decays [5j D + — > n~e + e + and D + — > K~e + e + . (Charge-conjugate 
modes are implicit throughout this Letter.) Short-distance FCNC processes in charm decays 
are much more highly suppressed by the Glashow-Iliopoulos-Maiani mechanism [f| than the 
corresponding down-type quark decays because of the range of masses of the up-type quarks. 
Observation of D + FCNC decays could therefore provide indication of non-SM physics or of 
unexpectedly large rates for long-distance SM processes like D + — ► n + V, V — > e + e~, with 
a real or virtual vector meson V. The LNV decays D + — > 7r~e + e + and D + — > K~e + e + are 
forbidden in the SM. They could be induced by a Majorana neutrino, but with a branching 
fraction only of order 10~ 30 . Any observation at experimentally accessible levels would be 
clear evidence of physics beyond the SM. Past searches have set upper limits for the four 
dielectron decay modes in our study that are of order 10~ 4 Q. The limits for corresponding 
dimuon modes are about an order of magnitude more stringent. 

The CLEO-c detector [H, 0, [l(l U| was used to collect a sample of 1.8 million e + e" 



^(3770) — > DD events (1.6 million D mesons) from an integrated luminosity of 281 pb _1 
provided by the Cornell Electron Storage Ring (CESR). From the interaction point out, 
CLEO-c consists of a six-layer low-mass drift chamber, a 47-layer central drift chamber, a 
ring-imaging Cherenkov detector (RICH), and a cesium iodide electromagnetic calorimeter, 
all operating inside a 1.0- T magnetic field provided by a superconducting solenoidal magnet. 
The detector provides acceptance of 93% of the full 4tt solid angle for both charged particles 
and photons. Charged particle identification (PID) is based on information from the RICH 
detector, the specific ionization (dE/dx) measured by the drift chamber, and the ratio of 
electromagnetic shower energy to track momentum (E/p). Background processes and the 
efficiency of signal-event selection are estimated with a GEANT-base d |l2l| Monte Carlo 
(MC) simulation program. Physics events are generated by EvtGen |l3( and final-state 
radiation (FSR) is modeled by the PHOTOS [1J| algorithm. Signal events are generated 
with a phase-space model as a first approximation of non-resonant FCNC and LNV decays. 

Candidate signal decays are reconstructed from well-measured charged-particle tracks 
that are consistent in three dimensions with production at the e + e~ collision point. Electrons 
with momenta of at least 200 MeV are identified with a likelihood ratio that combines E/p, 
dE/dx, and RICH information. Charged kaons and pions with momenta of 50 MeV or 
greater are selected based on dE/dx and RICH information. For each candidate decay of 
the form D + — > h ± e T e + , where h is either ir or K, we compute the energy difference AE = 
E ca ,nd — -^beam and the beam-constrained mass difference AM bc = \/ E^ cam — |p C and| 2 — Md+> 
where E C!mA and p C!m & are the measured energy and momentum of the h ± e T e + candidate, 
Ebea,m is the beam energy, and M D + is the nominal mass of the D + meson 0]. The resolution 
for these quantities is improved by recovering bremsstrahlung photons that are detected in 
the calorimeter within 100 mrad of electron trajectories. This provides a signal-efficiency 
increase of 13% — 18%, depending on decay mode. 

Events with D + candidates satisfying -30 MeV < AM bc < 30 MeV and -100 MeV < 
AE < 100 MeV are selected for further study. Within this region we define the "signal box" 
to be -5 MeV < AM bc < 5 MeV and -20 MeV < AE < 20 MeV, corresponding to ±3a 
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in each variable, as determined by MC simulation. The remainder of the candidate sample 
was used to assess backgrounds. 

The expected branching fraction for the long-distance decay D + — > 7r + — > 7r + e + e _ is 
within the sensitivity of this analysis (~ 1CT 6 ). We subdivide our candidates based on the 
mass squared of the final-state e + e~ (equal to the q 2 of the decay), with 0.9973 GeV 2 < 
m 2 +e _ < 1.0813 GeV 2 defining the 0-resonant region. We use this region both to veto the 
long-distance D + — > 07r + — > 7r + e + e~ contribution and to measure its branching fraction. 

Backgrounds in the D + — > h ± e T e + candidate sample arise from both DD and non-DD 
sources. In DD events double semileptonic decays are dominant. These typically have four 
or fewer tracks (including two real electrons) and large missing energy. Potential peaking 
backgrounds from three-body hadronic D + decays, such as K~tt + tt + , 7r + 7r + 7r~, and K + K®, 
are negligible because of the very small probability of misidentifying charged hadrons as 
electrons in CLEO-c (~ 0.1% per track), and because incorrect mass assignments result 
in AE outside the signal box. In non-DD events, including continuum e + e~ — > qq with 
q ^ c, r-pair events, and radiative return to ip(2S), non-peaking backgrounds arise from 
7-conversion and Dalitz decays of 7r° and r\. 

We have performed a "blind analysis." Signal-selection and background-suppression cri- 
teria were optimized using MC simulation before we open the signal box by minimizing the 
sensitivity variable 

, EZoC(n;N).V(n;N) 

where n is the observed number of events, N is the expected number of background events, 
C is the 90% confidence coefficient upper limit on the signal, V is the Poisson probability, 
N D + is the number of charged D mesons (as a function of integrated luminosity £), and e is 
the signal-selection efficiency The sensitivity variable S represents the average upper limit 
on the branching fraction that would be obtained from an ensemble of experiments if the 
true mean for the signal were zero. Sideband studies demonstrate that the MC simulation 
provides a good description of background events. 

Background associated with D semileptonic decays, mainly double semileptonic events 
with typically 4 or fewer tracks in the event with large missing energy (or semileptonic decay 
accompanied with 7-conversion or tt° and r\ Dalitz decays in the other side), is suppressed 
by a requirement on the energy Mother, the sum of the energies of all particles other than 
those making up the signal candidate. Small values of -Eothcr correspond to large values of 
missing energy in the event and are indicative of semileptonic decays in which neutrinos 
account for significant undetected energy. Optimization leads to different requirements on 
-Esther for different signal modes: -E thcr > 1-0 GeV for the 7r + e + e~ final state, -E ther > 1-3 
GeV for K + e + e~, and -E thcr > 0.5 GeV for the LNV modes if the number of tracks in the 
event is 4 or fewer. 

Background events from 7-conversion and from ir° and 77 Dalitz decays are suppressed by 
rejecting D + candidates with low effective dielectron mass. We use two kinds of dielectron 
effective mass squared variables for this purpose: m 2 +e - is computed for oppositely charged 
signal-side electrons and rh 2 e + e - is computed for all combinations of one signal electron with 
any unused oppositely charged track. We veto candidates if rn 2 +e _ < 0.01 GeV 2 or rh 2 + e - < 
0.0025 GeV 2 . 

The decay mode D + — > 7r + e + e" is susceptible to background from D + — > K Q s e + v e accom- 
panied by a semileptonic decay of the other D. This is suppressed by rejecting candidates 
when the signal n + and an oppositely charged track combine to give a mass M^^- that 
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FIG. 1: Scatter plots of AM^ C vs. AE obtained from data for each decay mode. The signal region, 
defined by -20 MeV < AE < 20 MeV and -5 MeV < AM hc < 5 MeV, is shown as a box. The 
two contours for each mode enclose regions determined with signal MC simulations to contain 50% 
and 85% of signal events, respectively. 



satisfies —5 MeV < M^+^- — M K o < 5 MeV, where M K o is the nominal Kg mass Q. 

After application of all background-suppression criteria, our intention was to eliminate 
multiple candidates (candidates in excess of one per mode per charge per event) by selecting 
the smallest |AM bc | among all that satisfy -5 MeV < AM bc < 5 MeV and -100 MeV < 
AE < 100 MeV. However, it turns out that there were no multiple candidate events. 

The residual background and the efficiencies after application of all selection criteria have 
been determined by MC simulation and are given for the four signal modes in Table |U The 
model used to describe FCNC and LNV decays is phase space. The efficiency is observed to 
be quite uniform over the Dalitz plot, with the exception of the two corners at low m e 2 e , which 
are depleted by the 200-MeV minimum-momentum requirement for electron identification. 

Scatter plots of AM^ C vs. AE for data events surviving all other cuts are shown in 
Fig. ^ For D + — > 7r + e + e~, two events lie in the signal box, with an expected background 
of 1.99. For all other FCNC or LNV modes there are zero events in the signal box. With 
no evidence of a signal, we calculate 90% confidence level (CL) upper limits (UL) on the 
branching fraction for each mode from the observed number of events (n) in the signal 
box, the signal-detection efficiency (e) , and the MC-estimated number of background events 
(N). We follow the Poisson procedure [7] to calculate the 90% confidence level coefficient 
(C(n; N)) upper limit on signal in the presence of expected background: 



UL 



C(n; N) 



(2-o-d 



+D~ 



(2) 



where o~d+d- [15] is the e + e~ — ► i/>(3770) — ► D + D~ cross section, £ is the integrated 
luminosity, and 2 • uo+d- • £ = 1.6 million is the number of charged D mesons in our data 
sample. Results are given in Table HJ we find no evidence of either FCNC or LNV decays. 
We separately measure the branching fraction for the resonant decay D + — > n + (f) — ► 7r + e + e~, 
finding two events in the signal region with an expected background of 0.04. 

Systematic uncertainties in these results can be divided into two categories: those related 
to background estimation and those arising from the signal-efficiency determination. 

For the background uncertainty, only the D + — ► 7r + e + e~ mode needs to be considered, 
as the other modes have zero observed events and the uncertainty in the expected number 
of background events does not affect their upper limits. For D + — > 7r + e + e _ , we compared 
the background estimate from the MC simulations with that of data from the AE-AM^ 
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TABLE I: Efficiencies (e), background estimates (N), observed yields (n), combined systematic 
uncertainties (<r sys t), and branching fraction results for four FCNC and LNV decay modes and for 
the resonant decay D + — > ir + <f) — ► n + e + e~ . Branching- fraction UL values are all at 90% CL. 



Mode 




N 


n 


O-syst (%) 


B (KT 6 ) 


n + e + e~ 


36.41 


1.99 


2 


8.7 


< 7.4 


7r~e + e + 


43.85 


0.48 





7.1 


< 3.6 


K+e+e- 


26.18 


1.47 





10.0 


< 6.2 


K~e + e + 


35.44 


0.50 





7.2 


< 4.5 


TT + (p(e + e~) 


46.22 


0.04 


2 


7.4 


2.7^ ±0.2 



sideband. The sideband estimate of the background in the signal box is about one standard 
deviation (a) higher than the MC estimate. Therefore our upper limit based on the MC back- 
ground estimate is conservative; the upper limit with the sideband-estimated background 
would be 5% lower. 

Sources of uncertainties that are common to all results are the number of D + (—3.2%, 
+4.5%), tracking (±1% per track or ±3% total), PID (±2.3%), and FSR (±4.0% for 7r ± e T e + , 
±3.3% for K + e + e~, ±3.5% for K~e+e + , and ±4.4% for -> TT+e+e", estimated by 
comparing the efficiency before and after bremsstrahlung recovery). 

Uncertainties in signal efficiency due to background-suppression cuts are estimated by 
comparing the efficiency before and after the cuts are applied: ±5.2% (7r + e + e~), ±1.1% 
(7r-e+e + ), ±7.3% (K+e + e~), ±1.0% (K~e+e+), and ±0.9% (tt+0 -> n + e + e-). 

Uncertainty from using the phase-space model (as a first approximation for non-resonant 
decays) for the FCNC and LNV signal efficiency estimation is assessed by (somewhat arbi- 
trarily) taking one quarter of the fraction of phase space which has non-uniform efficiency 
due to the electron identification momentum cut-off (200 MeV): ±2.8% (n ± e T e~) and ±3.8% 

For the results in Table HI we increase the upper limits to account for systematic uncer- 
tainties by decreasing the efficiency by lcr sys t (combined systematic uncertainty). 

In summary, we find no evidence for non-Standard-Model physics. There is no evidence 
either for the two rare (FCNC) decays or for the two forbidden (LNV) decays of charged D 
mesons to three-body final states with dielectrons. Finding no evidence for signals, we set 
90% confidence level upper limits: 

B(D + -> vr+e+e-) < 7.4 x 10~ 6 
B(D+ -> Tr~e+e+) < 3.6 x 10~ 6 
B(D + -> K + e + e-) < 6.2 x 10~ 6 
B(D+ -> K~e+e + ) < 4.5 x 10~ 6 

Our results for these dielectron modes are significantly more restrictive than previous lim- 
its, and reflect sensitivity comparable to the searches for dimuon modes 0]. Due to the 
dominance of long-distance effects in FCNC modes, we separately measure the branch- 
ing fraction of the resonant decay D + — > 7i + (f) — > 7T + e + e~, obtaining B(D + — > 07r + — > 
7r + e + e~) = (2.7l^*g ± 0.2) x 10 -6 . This is consistent with the product of known world aver- 
age 0] branching fractions, B(D + — > </>7r + — ► 7r + e + e _ ) = B(D + — > (f)iT + ) x B(<f) — > e + e _ ) = 
[(6.2 ± 0.6) x lO" 3 ] x [(2.98 ± 0.04) x 10" 4 ] = (1.9 ± 0.2) x 10~ 6 . 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
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